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The electronic contribution to the magnetically induced polarization in orthorhombic TbMnOs 
is studied from first principles. We compare the cases in which the spin cycloid, which induces 
the electric polarization via the spin-orbit interaction, is in either the b-c or a-b plane. We find 
that the electronic contribution is negligible in the first case, but much larger, and comparable 
to the lattice-mediated contribution, in the second case. However, we show that this behavior is 
an artifact of the particular pattern of octahedral rotations characterizing the structurally relaxed 
Pbnm crystal structure. To do so, we explore how the electronic contribution varies for a structural 
model of rigidly rotated MnOe octahedra, and demonstrate that it can vary over a wide range, 
comparable with the lattice-mediated contribution, for both b-c and a-b spirals. We introduce a 
phenomenological model that is capable of describing this behavior in terms of sums of symmetry- 
constrained contributions arising from the displacements of oxygen atoms from the centers of the 
Mn-Mn bonds. 



PACS numbers: 75.80.+q,77.80.-c 



I. INTRODUCTION 



Multifcrroic materials have been the subject of much 
excitement because they exhibit many interesting prop- 
erties and phenomenaiii£i£i^£ There are two basic sce- 
narios for the coexistence of two order parameters in a 
single phase. Either the two instabilities can both be 
present independently, or one can induce the other via 
some coupling mechanism. We are concerned here with 
magnetoelectric (ME) materials of the latter type, in 
which the primary instability is magnetic, and the re- 
sulting magnetic order induces an electric polarization. 
Such magnetically-induced (improper) ferroelectrics can 
be expected to display strong ME couplings, e.g., a strong 
dependence of the electric polarization on applied mag- 
netic field, or of the magnetization on the applied electric 
fields They could be extremely useful in many techno- 
logical applications, but most of the materials discovered 
to date either have too small of a ME coupling, or only 
operate at impractically low temperatures. Thus, it is 
essential to understand their coupling mechanisms more 
fully in order to design new materials with enhanced ME 
effects that can operate at higher temperatures. 

Among the best studied of the magnetically induced 
ferroelectrics are orthorhombic rare-earth manganites, 
0-i?Mn0 3 i^2i2i2iiaiiiiiSii2iii Intensive experimental and 
theoretical studies have clarified many questions regard- 
ing the origin of fcrroelectricity in these compounds. 
In HoMnC>3 and YMnC>3, the collinear E'-type anti- 
ferromagnetic (AFM) spin order induces a polariza- 
tion through the exchange striction mechanism ,L£i2 In 
TbMnC>3 and DyMnC>3, in contrast, the polarization ap- 
pears with the onset of spiral magnetic orde r 10 ' 11 as a 
consequence of the spin-orbit interaction.— 

On a microscopic level, the appearance of the po- 
larization can result either from a change in electron 
charge density that would occur even with ionic coor- 



dinates clamped (the purely electronic contribution), or 
from displacements of the ions away from their cen- 
trosymmetric positions as a result of magnetically in- 
duced forces (lattice contributions). In any theoreti- 
cal analysis of such materials, it is important to dis- 
tinguish between these two contributions and to calcu- 
late them separately in order to understand which mi- 
croscopic mechanisms are responsible for the appearance 
of the polarization. 8 ' 9 ' 12 ' 13 ' 14 Regarding the 0-RM11O3 
materials having the cycloidal spin structure, Xiang et 
al*2. and w o 12 ' 13 have demonstrated that in TbMnC>3 
the electronic contribution to the polarization is much 
smaller than the lattice contribution when the spin spi- 
ral is lying in the b-c plane. For this reason, our previ- 
ous wor k 13 ' 14 focused on a detailed analysis of the lattice 
contribution, while the electronic contribution was not 
studied carefully. The mode decomposition of the lattice 
contribution^ 3 - and its dependence on the spin spiral wave 
vector— revealed that the next-nearest-neighbor spin in- 
teractions are important not only for the formation of the 
spin spiral itself, but also for the induced polarization. 

Although in this particular case the lattice contribu- 
tion is dominant, it is not clear how general this result 
is. Picozzi et al£- showed that for orthorhombic HoMn03, 
in which the polarization is induced by collinear E-type 
AFM order, the electronic contribution to the polariza- 
tion is of the same order as the lattice contribution. 
Although the mechanism of polarization induction in 
HoMnC>3 is different from that in TbMnC>3, there is no a 
priori reason why the electronic contribution should be 
negligible in TbMnC>3. In fact, the first-principles study 
by Xiang et alH found that if the spin spiral lies in the 
a-b plane, the electronic contribution to the polarization 
is of the same order of magnitude as the lattice contri- 
bution. 

In this work, we study the polarization induced by the 
a-&-planc spin spiral and show that this case differs sig- 
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nificantly from the case of the 6-c-plane spiral. We focus 
mainly on the electronic contribution, analyzing it care- 
fully for both cases by considering how it varies for a 
structural model of rigidly rotated MnC>6 octahedra. We 
confirm the finding of Xiang et al*£ that the purely elec- 
tronic contribution to the polarization is not negligible 
for the case of the a-b spiral. Indeed, we find it to be 
quite sensitive to the choice of the calculation parame- 
ters, as well as on the octahedral rotation angles. Even 
for the case of the b-c spiral, we find that the electronic 
contribution can be quite significant if the octahedral ro- 
tation angles are varied away from the equilibrium values. 

We also construct a phenomenological model based on 
a symmetry analysis of the spin-orbit induced electronic 
dipoles associated with centrosymmetry-compatible oxy- 
gen displacements relative to the centers of the Mn-Mn 
bonds, finding that it is essential to take the Jahn- Teller 
orbital ordering into account from the outset. This model 
shows that b-c and a-b spin spirals need not be similar 
in terms of how the polarization is induced. Our work 
implies that the electronic contribution to the polariza- 
tion is generically expected to be much larger than was 
found for the specific case of the relaxed b-c spiral state 
in TbMn03, emphasizing the importance of considering 
both electronic and lattice mechanisms in any future the- 
oretical studies of this class of materials. 

The rest of the paper is organized as follows. In SccHTl 
we compare the electric polarization computed for fully 
relaxed TbMn03 with the spin spiral lying in the b-c 
and a-b planes for several values of the on-site Coulomb 
energy U in the LDA+U framework. In Sec. IIIII we fo- 
cus on the study of the purely electronic contribution to 
the polarization in the context of a structural model of 
TbMnC>3 in which the Mn and O ions form rigid octa- 
hedra. The dependence of the Berry-phase polarization 
on octahedral rotations is studied in Sec. IIII A[ while in 
Sec. IIII Bl we develop a symmetry-based phenomenologi- 
cal model in an attempt to explain the observed results. 
We discuss our findings in Sec. IIV[ and give a brief sum- 
mary in Sec. |Vj 



II. SPIN SPIRALS IN THE b-c AND a-b PLANES 

A spiral (or, more precisely, "cycloidal"— ) spin struc- 
ture forms in TbMnOs in the 6-c-plane below ~27K with 
the polarization lying along the c axis^ However, a suf- 
ficiently strong magnetic field applied along the b direc- 
tion causes the polarization to change its direction from 
c to a ("electric polarization flop"). It was suggested, 
and recently confirmed^ that this polarization flop re- 
sults from the change of the spin spiral from the b-c to 
the a-b plane ("spin flop"); the polarization simply fol- 
lows the spin spiral. We shall refer to these two magnetic 
states as the '&— c spiral' and l a—b spiral'. The former is 
incommensurate with a wavevector k s ~ 0.28, while the 
latter is commensurate with fc s = 1/4. 

In this section we first review the main results of our 



previous calculation s 13 ' of the polarization induced by 
the b-c spiral. We then present our new calculations for 
the case of the a-b spiral, and compare these two cases. 

We use a 60-atom supercell consisting of three Pbnm 
unit cells, corresponding to a spin-spiral wave vector of 
k s = 1/3, for both the a-b and b-c spirals. Although 
fc s = 1/4 experimentally for the a-b spiral, the use of 
four unit cells would be more computationally demand- 
ing, and the consistent use of k s = 1/3 facilitates com- 
parisons between the two cases. (An additional reason 
why k s = 1/3 is more convenient will be mentioned at 
the end of Sec. lUJBl ) 

Our electronic-structure calculations are carried out 
using a projector augmented- wav o 17 ' 18 method imple- 
mented in the VASP code packaged Since the local- 
density approximation (LDA) gives a metallic state for 
TbMnOs, we use on-site Coulomb corrections (LDA+U) 
in a rotationally invariant formulation^ The electric po- 
larization is computed using the Berry-phase method. 21 

In our previous work on the b-c spiral ) 13 ' 14 the struc- 
tural relaxation was performed in the absence of spin- 
orbit interaction (SOI), and we confirmed that no polar- 
ization is induced by the magnetic order in this case. In 
the presence of SOI, an electric polarization was found 
to appear. We decomposed it into electronic and lattice 
contributions by first keeping the ionic positions frozen 
at their centrosymmctric values while computing P, and 
then by repeating the calculation after allowing ions to 
relax. We found the electronic and lattice contributions 
to be -P G iec = 32^C/m 2 and Pi a tt = — 499/zC/m 2 respec- 
tively This result demonstrated that the lattice mecha- 
nism dominates over the purely electronic one for the b-c 
spiral in TbMn0 3 . 

The polarization values quoted above were calculated 
using U = 1 eV to match the experimental band gap^ 
We also studied the effect of the choice of U parameter 
(in a reasonable range of values from 1 eV to 4 eV ) on 
the induced polarization and coupling mechanism^ We 
found that while the absolute value of P becomes some- 
what smaller with larger U, the qualitative mechanism 
of polarization induction remains the same. 

In the case of the a-b spiral, we have now performed 
similar calculations of the polarization for the same set of 
U values from 1 eV to 4 eV. Table Q] shows the results for 
U = 1 eV and U = 4 eV. For these values we proceeded as 
in our previous work, taking a reference crystal structure 
that was fully relaxed in the absence of SOI, comput- 
ing the SOI-induced electric polarization p elec and ionic 
forces, and then using the latter, together with the com- 
puted Born charges and force-constant matrix elements, 
to predict P latt . (For intermediate U we only computed 
P elec , and did so in a simplified manner by using the ref- 
erence crystal structure that was relaxed at U = 1 eV, 
finding P olcc = 691 and 397/xC/m 2 for U = 2 and 3cV, 
respectively. These values are intermediate between the 
values for U = 1 and U = 4eV as expected.) 

For comparison, we also show in Table [J the results of 
similar calculations by Xiang et al.^ who used U = 2 eV 
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TABLE I: Purely electronic, lattice, and total polarizations 
along the c and a axes for the b-c and a-b spirals respectively. 
Results at Umu = 2eV (and also using Uxb = 6eV) from 
Ref. Q3 are shown for comparison. 
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-790 


740 
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-197 
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2 a 


331 a 


-462" 


-131° 



"From Ref. [H. 



on the Mn sites. However, their results are not directly 
comparable with ours, as they also included Tb / elec- 
trons, with /7xb = 6 eV on the Tb sites. In the b-c spiral 
case, as we mentioned before, the results do not depend 
strongly on the choice of U . Comparing in this case the 
results of Xiang et al. with our results for U = 1 eV, we 
find an agreement in that the purely electronic contribu- 
tion is negligible, and the total polarization values agree 
with each other within 10%. However, in the a-b spi- 
ral case there is no such agreement, which is perhaps not 
surprising in view of the very strong sensitivity of the po- 
larization to the value of U, as can be seen clearly in the 
table. It may also result in part from other factors, such 
as the different treatment of / electrons in the two calcu- 
lations, or the fact that they used a generalized-gradient 
(GGA) exchange-correlation while we used LDA. Nev- 
ertheless, a point in common is that both calculations 
predict that the electronic contribution is comparable or 
even larger than the lattice one for the a-b spiral. This 
leads us to conclude that the dominance of the lattice 
contribution that was found earlier for the case of b-c 
spiral is not a general phenomenon, but was special to 
that case. 

For the b-c spiral, the theoretical values in Ta- 
ble U are in satisfactory agreement with the value of 
~— 600/iC/m 2 found experimentally^ However, for the 
case of the a—b spiral, the comparison is more prob- 
lematic. Our computed polarization of 740/iC/m 2 
compares very poorly with the experimental value of 
~— 300/LtC/m 2 obtained by Yamasaki et al. for the re- 
lated Gdo.7Tbo.3Mn03 system^ However, as we shall 
discuss in Sec. IIV1 the polarization depends sensitively 
on the octahedral tilting angles, which may differ sig- 
nificantly for Gdo.7Tbo.3Mn03. Experiments on the a—b 
spiral in the TbMn03 system itself are somewhat am- 
biguous regarding both the sign and the magnitude of 
the polarization.— 

Note that the magnitude of the electronic contribution 
in the case of the a—b spiral falls rapidly with increasing 
U. Our calculations show that the band gap increases 
almost linearly with U. Fig. Q] shows the electronic con- 
tribution to the polarization for both a-b and b-c spi- 
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Average direct band gap (eV) 

FIG. 1: (Color online.) Dependence of electronic contribution 
to the polarization on the average direct band gap for the b-c 
spiral (circles, scale at left) and a—b spiral (squares, scale at 
right) when varying U. 



rals plotted versus the average direct band gap. One 
can see from the plot that the polarization is roughly in- 
versely proportional to the gap, up to a constant shift. 
A heuristic rationalization of this behavior can be given 
as follows. If we consider the derivative of the polar- 
ization with respect to ionic displacements, which is the 
Born effective charge, we know that this quantity can be 
expressed within density-functional perturbation theory 
in a Kubo-Grecnwood form involving a sum over terms 
that are inversely proportional to the differences of the 
eigenenergies of the unoccupied and occupied states^ 
The largest contributions are expected to come from the 
smallest energy denominators associated with states near 
the valence and conduction edges, so the overall sum 
should roughly scale inversely with the direct band gap. 
The same applies to other derivatives of the polarization, 
such as the dielectric susceptibility. If the derivatives of 
P have this behavior, it is not very surprising to find that 
the polarization itself has a similar behavior. 

In view of the results discussed above, the central ques- 
tion arises: Why are the cases of the a-b and b-c spiral 
so different? In the remainder of this paper, we attempt 
to shed some light on this question. For this purpose, 
we limit ourselves to a discussion of the purely electronic 
contribution to the polarization. We shall discuss at some 
length the dependence of the electronic polarization on 
atomic displacements, but only for displacement patterns 
that preserve inversion symmetry, such as those resulting 
from octahedral tilting in the Pbnm crystal structure. 



III. MODELING OF P clcc 

A. Structural model with rigid MnOs octahedral 
rotations 

To obtain a better understanding of the mechanism of 
the electronic contribution to the polarization, we con- 
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FIG. 2: (Color online.) Two initial configurations considered 
in the model of rigid MnOe octahedra. Rotations about c and 
b axes are indicated by white and dark blue curved arrows re- 
spectively, (a) Structure 1, which matches the physical Pbnm 
structure fairly closely, (b) Structure 2, with a fictitious pat- 
tern of octahedral rotations around the c axis only. 



sider a simplified structural model in which the crys- 
tal structure is composed of rigid corner-linked MnC>6 
octahedra, with the Tb ions remaining at their high- 
symmetry (0,0,1/4) Wyckoff coordinates. We then rotate 
the MnOg octahedra and study how the polarization de- 
pends on the rotation angles. All calculations within this 
model are done with U = 1 cV. While the values of the 
polarization computed with this U may not be realistic 
for the case of the a-b spiral, as discussed above, at this 
point we are interested in understanding the origins and 
behavior of the polarization, rather than making direct 
comparisons to experiment. Also, we want to compare 
the results with previous calculations ) 13 ! 14 most of which 
were done at U = 1 eV. 

Actually, before we even apply the rotations, we must 
first apply a Jahn- Teller (JT) distortion. Mn 3+ has a d A 
configuration in which the three majority-spin t2 g states 
are filled and the majority-spin e g levels are half-filled. 
The system is thus metallic in the absence of the JT dis- 
tortion; introducing it splits the e g levels and opens a gap, 
driving the system insulating. In our model, we take the 
JT distortion into account by pre-deforming the MnC>6 
octahedra such that the ratio of longest to intermedi- 
ate (along c) to shortest bonds lengths is 1.124: 1.004: 1, 
where these ratios have been extracted from our earlier 
first-principles calculations carried out with U = 1 eVi^ 3 - 

We then apply rotations to the octahedra. In Pbnm 
symmetry the rotations can be described as (a~a~b + ) in 
the Glazer notation^ meaning that out-of-phase and bi- 
phase alternating rotations occur around [110] and [001] 
axes respectively in the original cubic-pcrovskite Carte- 
sian frame. In the conventional framed used here, these 
correspond to the b and c axes respectively, and the spin- 
spiral wave vector propagates along b. 

In general, a different order of application of rota- 
tions leads to different final configurations, because ro- 
tations do not commute. Therefore, when describing the 
TbMnOa system in terms of MnOg octahedral rotations, 
one should carefully specify the meaning of the rotations 
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FIG. 3: (Color online) For b-c spiral, the dependence of the 
electronic contribution to the polarization on 6, the angle of 
rotation of the Mn06 octahedra about the b or c axis for 
Structure 1 or 2 respectively. (For the former, the rotation 
angle around c is 0/1.64; see text for details.) Symbols: first- 
principles calculations. Curves: result of the fit to the phe- 
nomenological model of Sec. IIII Bl 



and their order. For example, if we start with the ideal 
perovskite configuration and induce the Jahn- Teller dis- 
tortion, we can arrive at several possible initial configura- 
tions as shown in Fig. O If we apply (a~a~b + ) rotations 
to the configuration shown in Fig. [2fb) , regardless of the 
order of rotations, the final structure will not have Pbnm 
symmetry. However, applying a rotation around b fol- 
lowed by a rotation around c to the configuration shown 
in Fig.[^a), we will preserve the Pbnm symmetry. 

Fitting the angles of rotation to the relaxed structure, 
we find the rotation angles around b and c to be approx- 
imately 19.0° and 11.6° respectively. We thus constrain 
the ratio between these two angles to be 1.64, and treat 
the angle 8 around the b axis as the independent vari- 
able. We calculate the polarization for both the a-b and 
b-c spirals for a range of rotation angles —15° < 8 < 20°. 
We also made several calculations for the initial config- 
uration shown in Fig. [DJb), where the octahedra were 
rotated only around c. To distinguish between the two 
sets of calculations, we will refer to them as 'Structure 1' 
and 'Structure 2' respectively. 

The results of the calculations are presented in Figs. [3] 
and 01 Recall that all structures considered here have 
inversion symmetry, so that the Berry-phase calculations 
give us the purely electronic contribution to the polariza- 
tion induced by the SOI. These calculations reveal that 
even in the case of the b-c spiral, the electronic contri- 
bution to the polarization spans a wide range of values 
(^300/xC/m 2 ), depending on the octahedral rotations. 
This is yet another indication that this contribution is 
negligible in the relaxed b-c spiral structure only by co- 
incidence. 
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TABLE II: Classification of several quantities by their behav- 
ior under the mirror symmetries M x and M y : products of the 
spin components of two neighboring spins Si and S2, compo- 
nents of oxygen displacement vectors u, and components of 
the polarization vector P. 

Mx My 

+ 1 +1 Sl x S2x, Sl y S2y, SlzS2z U z U x , Uy, U 2 Z P z 

+ 1 —1 SlyS2z, Sl Z S2y Uy UyU Z Py 

— 1 +1 Sl x S2z, S\ Z S2 X U x U X U Z P x 

— 1 —1 Sl x S2y, Sl y S2 x U x U y 



FIG. 4: (Color online) For a-b spiral, the dependence of the 
electronic contribution to the polarization on 6, the angle of 
rotation of the MnC>6 octahedra about the b or c axis for 
Structure 1 or 2 respectively. (For the former, the rotation 
angle around c is #/1.64; see text for details.) Symbols: first- 
principles calculations. Curves: result of the fit to the phe- 
nomenological model of Sec. IIII Bl 



B. Phenomenological model 

To find out whether the observed dependence of the 
polarization on octahedral rotations can be explained 
within some relatively simple model, we decided to ana- 
lyze the possible contributions coming from each nearest- 
neighbor Mn-O-Mn triplet. Our notation is as follows. 
We use a, b, and c for the unit vectors in the global 
Cartesian frame of Fig. [2j We also attach a local Carte- 
sian frame to each Mn-O-Mn triplet as illustrated in 
Fig. [Ua), reserving z = 612 (the unit vector pointing 
from Mnl to Mn2), while x and y are chosen to form 
a right-handed triad with z such that x lies in the a-b 
plane. The origin of this frame is located in the middle of 
the Mn-Mn bond. The angle between &12 and the spin- 
spiral wavevector direction b is denoted by a, so that 
cos a = z • b. For the vertical bonds parallel to c, the 
local and global Cartesian frames coincide and a = tt/2. 
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FIG. 5: (a) Local (x,y,z) and global (a,b,c) coordinate frames, 
(b) Orbital ordering in TbMnOs. The d 3x 2_ r 2/d 3y 2_ r 2 or- 
bitals are aligned along the longest Mn-O bonds. Orbital 
order is uniform along the c axis. 



Our goal is to find dipole moments allowed by symme- 
try for each Mn-O-Mn bond, viewed in isolation. To find 
the total polarization, we need to transform these dipole 
moments back into the global Cartesian frame, average 
them over the spin-spiral period, and sum the contribu- 
tions from all Mn-O-Mn bonds in the 20-atom unit cell. 
One can show that there is no contribution to the polar- 
ization coming from the vertical bonds by using the fact 
that the magnetic moments on the Mn sites are collinear 
for these bonds. Therefore, we focus henceforth only on 
the bonds lying in the a-b plane. 

We expand the dipole moment of each Mn-O-Mn 
triplet in (i) bilinear products of spin components on 
the two Mn sites, and (ii) powers of oxygen displace- 
ments related to the M11O6 octahedral rotations. We 
invoke symmetry to determine the appropriate terms in 
this expansion, as follows. Consider the configuration 
shown in Fig.[^b). The Mn-O-Mn bonds have two mir- 
ror symmetries, M x and M y . Note that the bond does 
not have inversion symmetry because the JT distortion 
leads to an orbital ordering pattern, shown schematically 
in Fig. [5Jb) , which breaks M z . We emphasize that it is 
essential to take this JT distortion, and the associated or- 
bital order, into account. If instead one attempts to use a 
JT-frec pcrovskite structure as the reference for such an 
expansion, one finds the reference system to be metallic, 
so that the electric polarization cannot even be defined. 
Therefore, it is first necessary to establish an orbitally- 
ordered insulating state, and only then expand the elec- 
tronic dipoles in lattice displacements away from that 
state. Indeed, we initially attempted to derive a model 
built on the erroneous assumption of inversion symmetry 
for Mn-O-Mn bonds in a reference structure without JT 
distortions, but we could not arrive at a satisfactory low- 
order expansion that could simultaneously fit the data 
for both a-b and b-c spirals. 

We classify the products of spin components by their 
behavior under the two mirror symmetries, as tabulated 
in Tabic HTl We can then systematically expand the po- 
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larization in powers of displacements (u Xl u y ,u z ) as 

P x = A x °}S lx S 2z + AfJS lz S 2x 

+ [A x ^S lx S 2x + A$S ly S 2y + A£>S lz S 2z ]u x 
+ [A x ^S lx S 2y + A y 1 JS ly S 2x }u y 
+ {A£s lx S 2z + A^S lz S 2x ]u z + ... (1) 

Here we show only the terms that appear at zero and 
first order in u because the expression rapidly becomes 
tedious at higher order, but our analysis also includes all 
second-order terms. Similar expressions can be written 
for P y and P z . Projecting these contributions on the a, 
b, and c axes and averaging over the spin-spiral period 
and all Mn-Mn bonds in the unit cell, we arrive at 

P c (b " c) = sin(/){Co + C x u x + C z u z + C xx u\ + C yv u 2 y 

+ C zz u 2 z + C xz u x u z } (2) 

for the polarization in the case of the b-c spiral, and 

P^ a ~ b) = sin0{A o + A x u x + A z u z + A xx u 2 x + A yy u 2 y 

+ A zz u 2 z + A xz u x u z } (3) 

in the case of a-b spiral. The coefficients in Eqs. 

are just certain linear combinations of those appearing in 

Eq. (fT]) and the corresponding equations for P y and P z . 

The resulting expressions in Eqs. ([2][3]) for the polar- 
ization may be viewed as simple Taylor expansions in the 
oxygen displacements from the centers of the Mn-O-Mn 
bonds. However, some terms are missing because they 
arc forbidden by symmetry For example, terms linear in 
u y vanish after averaging along c because the contribu- 
tion from any bond in Fig. E^a) is canceled by the one 
from the bond above or below. The symmetry also im- 

plies that Pt c) = Pt C) = Pt b) = p c a ~ b) = after 
averaging over the spin-spiral period. The factor of sine/) 
comes from averaging the products of spin components 
over the spin-spiral period. In this model we did not 
consider next-nearest-neighbor spin interactions, which 
were shown to play an important role in the dependence 
of P c on the magnitude of the wave vector in the b-c spi- 
ral case i 13 ' 14 However, we argue that the contributions to 
the polarization coming from these interactions will have 
essentially the same form as Eqs. (J2JE]), but with sine/) 
replaced by sin(2(/>). In our particular case 4> = Tf/Z, so 
that sin</) = sin(2</)) and the inclusion of the next-nearest- 
neighbor interactions will only lead to a renormalization 
of the coefficients in the expansions. 

Eqs. ([2][3]) have no coefficients in common, showing 
that within this model there is no connection between the 
polarization in the b-c and a-b spirals. The results for 
each spiral case can be fitted independently with seven 
parameters, whose fitted values are given in Table IIIII 
In each of Figs. [3] and dj the resulting fits are shown as 
the solid and dashed curves that refer, respectively, to 
Structures 1 and 2 of Fig. O^a) and Fig.^b). Eqs. (TJED 



TABLE III: Fitted parameters d and Ai for the model of 
Eqs. (J2JI3J) for b-c and a-b spirals respectively. Units are 
/iC/m 2 , /iC/m 2 A, and /iC/m 2 A 2 for terms of overall order 
0, 1, and 2, respectively. Column labels are subscripts i. 








X 


z 


XX 


yy 


zz 


xz 




-285 


-3466 


211 


1297 


1143 


3263 


-30576 


Ai 


92 


-8245 


529 


555 


-1403 


361 


-33707 



clearly provide enough freedom to allow a very good si- 
multaneous fit to the results computed directly from first 
principles for both structures. If we go back and take the 
relaxed structures that were obtained directly from the 
first-principles calculations (for U=l eV and no SOI) and 
use the present model to evaluate the electronic contri- 
bution to the polarization, we obtain P c = 1600 fiC/m 2 
and P a = 31^tC/m 2 for the b-c and a-b spirals, to be 
compared with values of 1530 and 32/iC/m 2 computed 
directly from first principles, respectively. 

IV. DISCUSSION 

Our Berry-phase calculations of the electronic contri- 
bution to the polarization for various model TbMn03 
structures, in which rigid MnOg octahedra were rotated, 
show that the polarization in the a-6-spiral case behaves 
quite differently than for the fo-c-spiral case. Not only is 
the range of values different, but the qualitative depen- 
dence of P c icc on rotation angles is dissimilar. For the b c 
spiral, the polarization shows a parabolic dependence on 
the rotation angles, while for the a-b case the dependence 
is almost linear over a wide range of rotation angles. The 
phenomenological model considered above suggests that, 
at least from the point of view of symmetry, there is no 
relation between the b-c spiral and a-b spirals, so that 
the observed differences should not be very surprising. 

Considering that the octahedral rotation angles and 
oxygen displacements become quite large, we obtain quite 
good fits of the dependence of the polarization on rota- 
tion angles from our expansion of the dipolcs on oxygen 
displacements away from the Mn-Mn bond centers. The 
observed behavior results from the fact that the coeffi- 
cients in front of u x and u z are much larger in the a-b 
spiral, while the quadratic coefficients in front of u xx and 
u zz are much larger in the b-c spiral (see Table IIII[) . It 
still remains to understand how the most important co- 
efficients in the expansion get their values based on some 
microscopic model of bond hybridization. 

If wc compare our computed polarization for the 
relaxed TbMn03 structure in the a-b spiral case 
(740/uC/m 2 , see Table HJ to the experimental value for 

Gdo.rTbo.sMnOs in Ref. [H ( 300/zC/m 2 ), wc find 

very poor agreement. However, since Gd has a larger 
radius than Tb, the MnOe octahedra will be less tilted, 
reducing the electronic contribution to the polarization 
(see Strucure 1 in Fig. 3]). This effect may help explain 
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the observed difference. However, it should also be kept 
in mind that the strong sensitivity of the polarization 
to the choice of U in the case of the a-b spiral means 
that any prediction of the polarization made within the 
LDA+U framework will have a much larger uncertainty 
than for the b—c spiral case. The use of linear-response 
techniques to compute the effective parameters for the 
LDA+U method^ could thus be appropriate here. How- 
ever, the very use of the LDA+U method itself may be 
questionable, and it may be worth exploring the suit- 
ability of other methods, such as GW quasiparticle^ or 
dynamical mean-field theory^ approaches, for comput- 
ing the polarization in this case. 

We have seen that the octahedral rotations can sig- 
nificantly change the polarization. Although we have 
focused here only on the electronic contribution, it ap- 
pears likely that the lattice contribution will also depend 
strongly on rotation angles. Such a calculation of the 
lattice contribution is problematic, however, because one 
wants to consider the SOI-induced symmetry-breaking 
distortions away from a reference structure that is not 
itself an equilibrium structure in the absense of SOI. 
In principle it may be possible to compute these using 
an approach similar to that in Ref. Il3l . That is, one 
would compute the force-constant matrix and dynami- 
cal charges (in the absence of SOI) and the SOI-induced 
forces for a given configuration of octahedral rotations, 
and use these to predict the induced amplitudes of the 
infra-red-active phonon modes and the resulting lattice 
contribution to the polarization. We have not pursued 
such an approach here, as it would take us beyond the 
intended scope of the present work. 



V. SUMMARY 

We have used first-principles methods to compute the 
electronic and lattice contributions to the spin-orbit in- 
duced electric polarization in the cycloidal-spin com- 
pound TbMn03 with the spin spiral in the b c and a b 
planes. In the latter case we find that the electronic con- 
tribution is of the same order of magnitude as the lattice 
contribution, in strong contrast to previous studies of the 
b-c case. 

We have studied the electronic contribution to the po- 
larization in detail by considering a structural model 
based on rigid rotations of MnOg octahedra. We have 
shown that the electronic contribution to the polarization 
can change significantly with rotation angle even in the 
case of the b—c spiral, thus demonstrating that our previ- 
ous neglect of this contribution was justified only because 
of an accidental property of the relaxed Pbnm structure. 
We have introduced a phenomcnological model that ex- 
pands the electronic contribution to the polarization up 
to second order in the oxygen displacements from the 
Mn-Mn midbond positions, and have shown that it can 
explain the quite different behavior of the polarization in 
the b-c and a-b spiral cases. 
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